Effective Operator for dd Transitions in Nonresonant Inelastic X-ray Scattering 
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Recent experiments by Larson et al.^ demonstrate the feasibility of measuring local dd excita- 
tions using nonresonant inelastic X-ray scattering (IXS). We establish a general framework for the 
interpretation where the dd transitions created in the scattering process are expressed in effective 
one-particle operators that follow a simple selection rule. The different operators can be selectively 
probed by employing their different dependence on the direction and magnitude of the transferred 
momentum. We use the operators to explain the presence of nodal directions and the nonresonant 
IXS in specific directions and planes. We demonstrate how nonresonant IXS can be used to extract 
valuable ground state information for orbiton excitations in manganite. 

PACS numbers: 78.70.Ck, 71.20.Be 



I. 



INTRODUCTION 



Transition-metal compounds display a wide variety of 
exciting phenomena such as high- Tc superconductivity, 
colossal magneto resistance, and metal-insulator tran- 
sitions resulting from the strong interplay between the 
charge, spin, and orbital degrees of freedom. Knowl- 
edge of the electronic and magnetic structure can be 
obtained by a variety of spectroscopic techniques, such 
as optical spectroscopy, photoemission. X-ray absorp- 
tion. A general understanding of spectroscopy and its 
uses has always been a crucial aspect in the advance- 
ment of condensed-matter physics. Recently, inelastic X- 
ray scattering (IXS), both on and off resonance, has at- 
tracted considerable attention. Resonant IXS (RIXS)^ is 
a second-order process dominated by p- A, where p is the 
momentum and A the vector potential. The incoming X 
ray excites an electron from a deep-lying core state into 
the valence shell, and one measures the radiative decay of 
the core hole. For RIXS at the transition- metal L and M 
edges, the dipole transitions create dd excitations that in- 
clude spin-flips due to the strong intermediate-state spin- 
orbit coupling^ ^. The transitions can be described with 
an effective operator approach^ in the fast-collision ap- 
proximation. RIXS at the K edge, where excitations are 
predominantly shake-up processes by the strong Is core 
potential^, can be related in certain limits to the dynamic 
structure factor Sc^{(jo)^ . 

Nonresonant IXS, on the other hand, involves the in- 
teraction 
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with q = k— where e/e' and k/k' are the polariza- 
tion and wave vectors of the incoming/outgoing X rays, 
respectively. The interaction due to the term is in 
principle weak. However, the successful experiments by 
Larson et al.^ demonstrate the feasibility of measuring 
local dd excitations and comparison of the radial ma- 



trix elements indicate that these transitions can in prin- 
ciple also be observed in other transition-metal and rare- 
earth compounds. The appealing feature of nonresonant 
IXS is that the cross section is proportional to the dy- 
namic structure factor: ^ Sc^{<jj)^ which through 
the fluctuation-dissipation theorem is connected to the 
imaginary part of the density response function. 
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with Eg the ground state energy and where the Hamilto- 
nian for transition-metal compounds, includes strong 
many-body interactions; pq is the density operator 
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where cj^^^ creates an electron in a state V^kncr where 
n is the band index and a =t,i- Without matrix ele- 
ments pq becomes equivalent to the charge density op- 
erator pq = Xlkcr ^k+q cr^kcr- Howcvcr, matrlx elements 
play a crucial role in the understanding of the inelastic 
scattering. Using Wannier functions, we can rewrite the 
matrix element in terms of scattering from an atom at 
site Ro to a site R, 



(V^k+q,nv|e*'*'''|V^kna) 
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where (p-Raa is a localized Wannier orbital of type a at 
site R. The operator pq can create charge-transfer tran- 
sitions (with R 7^ Rq), plasmon excitation^^HH^ ^nd 
dipolar and higher-order transitions from core to va- 
lence states (where R = Rq)"*^^ ^'^ - Recently, it was 
demonstrated^ that it is also possible to measure dipo- 
lar forbidden dd transitions. By tuning the transferred 
momentum^ one is able to maximize the intensity for 
local transitions (R^ = Rj) within the shell, for which 
the matrix element is 

(^k+q,n'k*'*'''l</^kn) = ((/>3d,a' 
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where ol denotes the diflFerent orbit als. This approx- 
imation assumes reasonably weh 'locahzed' Warnier or- 
bitals. These sharp dipole-forbidden transitions were first 
observed by Larson et al} for large q in NiO and CoO. 
The use of large wave vectors would allow a resolution 
of 30 meV or less, making nonresonant IXS a promising 
tool to study local crystal field and orbiton excitations. 
It could be used, for example, to study the Jahn- Teller 
distortions in manganites and the nature of the small 
crystal field distortions in early transition-metal oxides. 
Although the angular dependence in NiO and CoO has 
been analyzed with density-functional theor}/^ and small- 
cluster calculations^^, the detailed nature of the angular 
dependence and how to use it is not well understood. In 
this paper, we express pq as an effective operator, derive 
the selection rules governing the local dd transitions, and 
give explicit angular dependencies. We provide an expla- 
nation for the remarkable intensity variations in certain 
directions. The derivation for the angular distribution of 
the intensities is partially based on symmetry rules and 
should therefor be very generally valid. We show how IXS 
can be used to extract detailed ground-state information 
by treating the IXS from orbitons (excitonic orbital ex- 
citations not involving the Hubbard U) in manganites. 



TABLE I: The angular dependence t/mm'(q) = (m^|pq|m) of 
the scattering between orbitals m and m! . The real orbitals 
are Z^'* — VSxy, VSyz, VSxz , |>/3(^^— form — 

—2, —1, 0, 1,2; q = (x, y, z) — (sin^cos(/?, sin^sin(/?, cos^), in 
conventional spherical coordinates and Lp. For off-diagonal 
matrix elements (m ^ m') for t2g {m — ±1, —2) orbitals, m" 
denotes the t2g orbital for which m" / m^m . The coordinate 
Tm — y^x^z for m = 1, —1, ±2. 

m m Umm' (q) 



^2(-p^ + I) + |^4(35z^ - 30z^ + 3) 
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±1 [-1^2 + 1^4(7^2 -3)]2i'i' 

±1 2 [T^^2 + |V3/l4{7(x2-y2)±2}]z£> 

±2 [A2 + IAa{7z'^ -l)]z''il 
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and an angular dependence Z^^(q). Of the summation 
over L only the values (monopolar), 2 (quadrupolar) , 
and 4 (hexadecapolar) remain. The factor 



II. EFFECTIVE OPERATORS AND 
SELECTION RULES 

For local dd transitions, it is convenient to express the 
3(i Wannier functions (b-Rnm = 

and r = r/r, in terms of a radial function Rsdi'^) and a 
real angular part 
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where m = |m|, Nm 

0, (— l)^sgn(m) for m = and m > 0, respectively; 

is a spherical harmonic. The functions (y) 
are known as tesseral harmonics and convenient when 
dealing with transition-metal compouns since the val- 
ues m = —2,-1,0,1,2 correspond to the 3d orbitals 
dxyi dyz, dsz'^-r^^ d^z^ dx2_y2^ respectively, see Table[l| We 
can expand the exponent in Eqn. ([T]) in terms of Bessel 
functions and tesseral harmonics, 
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where M = — L,— L + 1,-- - ,L and Jl is a Bessel function 
of order L. Spherical-tensor algebra gives an effective 
transition operator 
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where C\^^^ 1 ^ are Clebsch-Gordan coefficients, and 

limi,l2Tn2 ' 

Pl{q) = / drr'^Rsd{^)jL{Q'^)R3d{'^) is the reduced matrix 
element of the Bessel function. For brevity, we implic- 
itly assume the dependence on q and q in the remainder. 
Al = Po, ^^2, 1^4, for L = 0, 2, and 4, respectively. In 
second quantization, the transition operator is 
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where d}^^ creates an electron with spin a =T, j in the 
3d orbital with index m. The transition probability is 



in 1)1 ■ ^m±,sgn(mM)|m±M| 
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with P^^ 2sm, 2sm for m ± M > and m ± M < 0, 
respectively, and = 2 - ^m,o^M,o. 

For the monopole term (L = 0), the scattering is elastic 
and isotropic, since 



^ ^ d\n(j^rna-, 
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and 4°^(q) = 1. 

For inelastic scattering in transition-metal systems, 
the coefficient contains the simple selection rule 



m! = sgn(mM)|m ± M|, 



(12) 



consisting of a one-particle transition operator 



M 



probed by q through a reduced matrix element AL^q) 



under the conditions that \m'\ < 2 and sgn(mM) = 1 
for m' = 0. This selection rule helps us to obtain an 
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understanding of nonresonant IXS. As an example, let 
us consider a Cu^+ ion in D^h symmetry for quadrupo- 
lar {L = 2) scattering. The ground state is \d^2_y2) 
(m = 2), where the underline indicates holes. When 
measuring along the [001] direction, the only nonzero an- 
gular term is Zq^^ = ^z'^ — ^ (^m = 0), where we use 
q = {x^y^z) = {s'mO cosip^s'mO simp^cosO) in conven- 
tional spherical coordinates and cp. There is no inelas- 
tic scattering, since the relevant transition operator Wq 
only contributes to the elastic intensity (m = 2 ^ 2). 
When measuring with q in the xy plane, Z^"^ is zero for 



with 



odd M. For even M, 



gives elastic scattering; wz 



does not contribute since transitions to = —0, —4 
are forbidden since sgn (mM) = 1 is not satisfied for 
m' = and \m'\ < 2 for 3d electrons. The only in- 



elastic scattering is due to W2 giving d^2_ 
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(or m = 2 ^ 0) with a zf^ = ^(x^ - f) angular de- 
pendence. In addition, we can easily see that transitions 
d^2_y2 d^y (m = 2 ^ —2) are possible for operators 
with symmetry —0 and —4 using the inverse relationship 
M — sgn (mm') I m ± m'\. However, —0 does not sat- 
isfy the condition that sgn(mm') = 1 for M = and is 
therefore not allowed. The transition d^2_„2 
therefore hexadecapolar (L = 4, M = —4). 

Although all the off-diagonal terms for the quadrupo- 
lar scattering are determined by a single M, often two 
different M values interfere. It is therefore convenient to 
define a total angular dependence between two 3d orbitals 
by f/mm'(q) = (^^Pql^), which are given in Table [l[ We 
now demonstrate how the selection rules and symmetry 
can help in extracting valuable ground-state information. 



III. 



IXS IN HIGH-SYMMETRY DIRECTIONS 



One of the striking features of nonresonant inelastic 
X-ray scattering in NiO is the presence of a [001] nodal 
direction absent in CoO. Larson et al} ascribed this to a 
"q-selection rule" associated with the nearly cubic point 
group symmetry of NiO. The absence of inelastic scatter- 
ing, (/|pq|^) = 0, implies Pq|^) = const |^). This occurs 
when H and pq commute, [H^ pq] = 0. In general, this 
rarely happens, but Pq can commute with parts of the 
Hamiltonian, which then might give rise to nodal direc- 
tions for all excited states if Pq|^) = const |^) and oth- 
erwise gives a nodal direction for the excited states for 
which (/|pq = (/I const. Below we will first discuss the 
commutation relations between pq and the crystal-field 
operator and then those with the Coulomb interaction. 
For q = gz, the only nonzero angular dependence Zq^^ 
selects the operators Wq . This gives a density operator 
diagonal in m. 
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In several common crystal- fields, such as Oh and D^h^ 
m is a good quantum number and therefore Pq(z) com- 
mutes with the crystal-field operator. Along [100], Pq^ 
is no longer diagonal but, since it is a unitary transfor- 
mation over 90° of p^z, still only contains e ^ e and 
t2 ^2, but no e ^ ^2 scattering. Eigenstates of the 
octahedral crystal field therefore scatter among them- 
selves. Thus, for a system with a Oh or 1)4^ symmetry 
a nodal direction can occur when q || C4, where C4 is 
a fourfold symmetry axis. Since in general pq does not 
commute with the Coulomb interaction, the eigenstate 
of the crystal field should also be an eigenstate of the 
Coulomb interaction. This occurs for high-spin 3d'^ con- 
figurations, except 3d^ and 3d^ . For example, the ground 
state for a Ni^+ ion (in the absence of spin-orbit cou- 



ground and final state orbitals: 



^^^^ ^^^^ ^^^^ ^^^^^ 

^>lfie!) lOifie!) i^^,>|fi4) i^+)|ti^'!> 

a^=li\a^=0 0.9||0.1 0.75i|o.25 0.5| 1 0.5 
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FIG. 1: (color online) Inelastic scattering from e ^ e* in the 
manganites. (a) The top part shows the lowest e orbitals for 
|e+) = a\2) + p\0) for = 1,0.9,0.75, and 0.5. Excitations 
are made into the orbital |e^) = a\0) — f3\2). The angular 
dependence for = 0,0.1,0.25,0.5 is the same as those for 
1 — a^. The inversion of the relative energy positions of the 
orbitals does not affect the nonresonant IXS, which is propor- 
tional to |(e±|pq|e±)|^. (b) The angular distribution of the 
nonresonant IXS intensity in the quadrupolar region {A2 = 1 
and A4 — 1 — A2 — 0). The intensity in a certain direction is 
proportional to the distance to the origin. The insets show the 
intensity in the xy plane. The intensities give the sum over 
both orientations of the e± orbitals; (c) idem for A2 = 0.67; 
(d) idem for A2 = 0. 
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pling) is d^z^_r2^d^2_y2^. This is an eigenstate of the 
crystal field, the Coulomb interaction, and pqz- There- 
fore, no inelastic scattering occurs along all six C4 axes. 
On the other hand, for Co^+, the ground state is given by 
1^) = a|^,^^d^^(4Ti)) +/3|4^4^(4Ti)), where the mix- 
ing occurs due the Coulomb interaction. This is not an 
eigenstate of the octahedral crystal field and the oper- 
ators pgz, and inelastic scattering in the [001] direction 
therefore occurs into the excited multiplet 



l/)=/3|d*,T^eT('ri)) 



(15) 



which is about 2.4 eV higher in energjP^. Also, devia- 
tions from a ^26^ ground state due to band effects and 
a lowering of the crystal fielcF^ give rise to inelastic scat- 
tering. Spin-orbit coupling can also remove a nodal di- 
rection, but its effect is small except when the spin-orbit 
coupling lifts a degeneracy. Finally, the conditions for 
nodal directions above are generally not satified for low- 
and intermediate-spin ground states. 



IV. ORBITON EXCITATIONS IN 
MANGANITES 

Orbital physics plays an important role in manganites 
which have been extensively studied for their magnetore- 



sistive behavioi^. LaMnOs is known to have an alter- 
nating dsx^_r2 /dsy2-r^ orbital ordering, and half-filled 
systems often display the CE type structure with charge 
and orbital order occurring in an unconventional zigzag 
magnetic structure. However, others have contested this 
ionic picture stating that the ground state is more com- 
plex. For example, for Lao.5Sri.5Mn04, it was claimed, 
based on X-ray magnetic linear dichroism experiments^, 
that a significant out-of-plane character was mixed in, 
giving an orbital ordering close to dx2_^2 /dy2_^2 . We 
demonstrate the extreme sensitivity of the angular de- 
pendence of nonresonant IXS on the detailed nature of 
the ground state. The ground state for Mn^+ is given 
by a d^^^de^ configuration. Degenerate e orbitals ^322-^2 
(m = 0) and dx2_y2 (m = 2) are sensitive to distortions 
leading to lowest states given by \e±) = ±a|2) + /3|0), 
with + = 1. The phases account for the different 
orientations of the orbitals due to orbital ordering often 
found in manganites. With nonresonant IXS, one can 
make excitations into the empty state |e5-) = a\0) ^ P\2) 
by exciting an electron from the e {d^^^de^^i d^^^de^^) 
or t2 states {d^^^de^^ d^^^de^^de^^). For the former, 
the inelastic intensity summed over both e± orientations 
is 



2a^/3\Uoo-U22f+2{a^-0'rUi 
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FIG. 2: (color online) Nonresonant inelastic scattering 
for the manganites with an electron in the |e) states as 
in Fig. Qa), but now for scattering t2 e* for = 
1,0.9,0.75,0.5,0.25,0.1 and 0. Note that now and 1 - 
are inequivalent. (a) Nonresonant IXS in the quadrupolar re- 
gion (A2 — 1 and Aa — 0); (b) idem in the hexadecapolar 
region {A2 = and A4 = 1). The insets show the scattering 
intensity in the xy plane. 



see also Table [l[ Figure [T] shows the angular dependencies 
for different ground-state orbitals given by a^, see Fig. 
[TJa). Let us first look at the quadrupolar region {A2 — 1, 
A4 = 0), see Fig. [lib). The angular dependence can be 



rewritten as Sa^/^^p^ _ i)2 + 2(^2 - ff^f^x^ - ff 
and can be straightforwardly used to extract the value of 
a. A typical experiment is comparable to those of Lar- 
son et al}. One fixes the magnitude of the transferred 
momentum and changes its angle. For example, the in- 
tensity when rotating the angle from the [001] to the [100] 
direction depends strikingly on a. For a = 1, which cor- 
responds to an electron in the — y'^ orbital, the IXS 
intensity is predominantly along the x and y directions. 
For ^ = ^1 the angular intensity is predominantly along 
the z direction. However, when the hexadecapolar con- 
tribution increases, such an interpretation is less straight- 
forward, see Fig. [TJc) and (d). However, simplifications 
occur when looking in specific directions or planes. First, 
note that the [001] direction is only nodal when |e+) = |2) 
or |0), which are good eigenfunctions of pqz. However, 
more quantitative information can be obtained by com- 
paring the intensities along [001] and [100] or [010]. The 
e ^ e scattering along [001] depends on the parameters 
Qo = -M + 6A4 and (32 = ^2 + ^4, see Eqns. ([l3| 
and (14). However, for inelastic scattering, we can re- 



move one of those parameters by rewriting the scattering 

operator as = Qo^e + Em7^0,a(Qm - Q^)d^mcj^ma - 

The number operator ng only gives elastic scattering, 
so effectively the e ^ e scattering only depends on 
Q20 = Q2 - Qo = 2A2 - 5A4. Along [100], only 
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contains e ^ e and t2 ^2, but no e ^ ^2 scatter- 
ing. The e ^ e scattering again only depends on Q2o- 
A straightforward calculation gives for the ratio of the 
nonresonant scattering intensities in the [001] and [100] 
direction 

/[ooi]//[ioo] = 16a2/3V(3 - Sa^p^). (16) 

Note that, since = 1 — a^, this ratio is independent of 
the reduced matrix elements and can be used to extract 
directly the nature of the Jahn- Teller distorted state. 
From the factor Q20 = 2^2 — 5^4, one also sees a destruc- 
tive interference independent of a between quadrupolar 
and hexadecapolar terms in the and z directions, 

which is clearly visible in Fig. [ijc), where the intensity 
in those direction is almost zero. A clear feature that 
displays the change in ground state as a function of 
is the direction of the dominant lobes in the xy plane, 
see insets in Fig. [TJb)-(d). Again, this can be used to 
obtain information on ground-state properties. The in- 
tensity in the x or y direction ^Q2q{3 — 8a^ p'^) decreases 
with increased mixing of the — y'^ and — or- 
bitals, whereas along the [110] directions, the intensity 
2a^/3^(A2 + ^A^)'^ increases. 

Although the measurement of the e ^ e* orbiton ex- 
citation has the great advantage that the energy loss can 
be directly related to the energy for an orbiton excita- 
tion. Since this energy is of the order of 0.3-1.5 eV, 
this feature might in practice be difficult to distinguish 
from the elastic line. However, valuable information can 
still be obtained from the study of the ^2 e* exci- 
tation, where an electron from the t^g spin is excited 
into the empty e* orbital. This excitation is at higher 
energy loss due to the crystal-field splitting between e 
and ^2, which is of the order of 2-2.5 eV. The angular 
dependencies are shown in Figure |2j First, note that 
the spectra for and \ — o? are no longer equivalent. 
Second, the x, and z directions are all nodal since 
no ^2 ^ e scattering occurs. Simplifications also oc- 
cur in the xy plane. Scattering from |0) |±1), | — 2) 



requires M = ±1,-2 and |2) 1^1), I — 2) requires 
M = ±1(±3) and (—4), respectively (M- values in paren- 
theses give pure hexadecapolar terms). All the odd M 
are zero in the xy plane and only |0) |— 2) (M = —2) 
and 1 2) |— 2) (M = —4) remain, giving an angular 
distribution in the xy plane 2a^(/7o,-2)^ + 2/5^ ([/2, -2)^. 
The relative strengths of /7o,-2(^, ^, 0) = ^Q2oxy and 
^2,-2(^,^,0) = ^A4xy{x'^ — y'^) refiect the amount of 
3^2 _ ^2 ^2 _ y2 character in the ground state, re- 
spectively. Note that in the quadrupolar region, we have 
a simple scaling of the xy dependence, see Fig. [2|a). The 
advantage of using hexadecapolar excitations, see Fig. 
|2|b), is the clear change in angular dependence from a 
four-lobed to an eight-lobed shape as a function of a. 
V. SUMMARY 

The nonresonant inelastic X-ray scattering for local 
dd transitions has been analyzed. The strong sensitiv- 
ity of the angular dependence on the detailed nature of 
the ground state in combination with the experimental 
degrees of freedom (scattering angle, incoming energy) 
and the possible high resolution make nonresonant IXS 
a powerful tool to study crystal-field and orbital excita- 
tions. Future theoretical work should include an analysis 
for the rare-earths. 
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